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Abstract: 

We calculate the complete 0(a s ) corrections to the quark decay b —> ccs taking full account 
of the quark masses, but neglecting penguin contributions. For a c to the b quark mass 
ratio m c /m & = 0.3 and a strange quark mass of 0.2 GeV, we find that the next-to-leading 
order (NLO) corrections increase F(b —> ccs ) by (32 ± 15)% with respect to the leading 
order expression, where the uncertainty is mostly due to scale- and scheme-dependences. 
Combining this result with the known NLO and non-perturbative corrections to other B 
meson decay channels we obtain an updated value for the semileptonic branching ratio of 
B mesons, B S l, of (12.0 ± 1.4)% using pole quark masses and (11.2 ± 1.7)% using running 
MS masses. 


Submitted to Phys. Lett. B. 





1. Owing to the newly developed tool of expanding in the inverse heavy quark mass [T|] , 
the theoretical description of weak inclusive decays of heavy mesons now rests on more solid 
ground than ever. Since in such decays the energy release is large compared to the masses of 
the final state particles, the process takes place essentially at small distances and, in leading 
order in the heavy quark expansion (HQE), is described by the underlying quark decay, 
ffadronic corrections only enter at second order in the ffQE and should be 1 GeV 1 2 /ml, 
which is around 5% for B decays with a b quark mass mb ~ 5 GeV. Thus the accuracy of 
theoretical predictions for hadronic quantities like, say, the semileptonic branching ratio is 
not so much limited by our necessarily incomplete knowledge of (non-perturbative) hadronic 
matrix elements, but rather controlled by our knowledge of perturbative corrections to the 
free quark decay. 

This issue has recently attracted much attention in connection with the summation of 
certain terms of the perturbative series, namely the asymptotically leading ones of order 
aT'PS 0- Although this program can straightforwardly be applied to semileptonic B 
decays^] [[j, there are severe problems (both technical and conceptual ones) with applying 
it to nonleptonic channels, so that in this letter we only deal with the first order radiative 
corrections. 

Until recently, full 0(a s ) corrections were only known for the semileptonic decay b —> ceis 
1, 0 and for b — > ctu [Q. Although it is known that the exchange of gluons between 
quarks of unequal masses can yield big effects (cf. the extreme case of an infinitely heavy 
quark investigated in Ref. 0), finite c quark mass effects in the 0{a s ) corrections to the 
nonleptonic decay b —> cud and a rough estimate^ of c quark mass effects in b —> ccs have 
only recently been obtained [|5|. In Ref. |TUJ part of the effects of finite charm and strange 
quark masses in the radiative corrections to b —» ccs were taken into account, based on the 
calculation done in Ref. [Tij]. In this letter, we complete the calculation of finite quark mass 
effects in the 0(a s ) corrections to b —> ccs, neglecting penguin corrections. We exploit this 
result to give an updated prediction for the semileptonic branching ratio Bsl of B mesons. 


2 . In calculating the decay rate T(6 — ► ccs), we start from its representation as the 
imaginary part of the relevant forward-scattering amplitude: 

T(6 —■> ccs ) = — Imi [d^x ( b\T £j) c=2 (x) £w C=2 (0) 16). (1) 

mb J 

£-w' the effective Lagrangian that describes the decay process in the limit of an infinite 
W boson mass and to first order in the weak coupling. It can be written as 

■C c * 2 = - % v c. Ki £ CiMOiM. (2) 

V £ i =i 


In this letter we conform to the notation of 0, where the operators 0\ and 0-2 denote 
current-current operators with a colour-non-singlet and colour-singlet structure, respec¬ 
tively, whereas the remaining operators are due to the admixture of penguin contributions. 


1 Cf. H for an explicit calculation of the a 2 /3o term in the decays b —> uev and b —> ceu. 

2 See also ||. 
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The Cj are perturbatively calculable short-distance coefficients that describe the physics be¬ 
tween the scale of the W boson and the characteristic hadronic scale of the process, which 
is of the order of the b quark mass. If we neglect the effects of strong interactions, then 
only Oo has a non-vanishing Wilson-coefficient. In next-to-leading order (NLO), both the 
coefficients <7 and the operators O t depend on the scheme used to deal with 75. To that 
accuracy, the decay rate T (6 —> ccs ) can be written as 

5 6 i 

T(b ccs ) = \V cb \ 2 \V cs \ 2 PA(x c ,x c ,x s ) (3) 

i=l j=1 


G 2 F ml 

647T 3 


| v cb \ 2 \ Vcs\ 1 PH(x c , x c , x s )k(x c , x s , h) K(x c , x s ,h), 


( 4 ) 


where dij = 1 + \a s (mw) — a s (^)]/ 7 r + Ay,- a s (fi)/7T. The function n is defined in such a 
way as to contain the LO effects, whereas the product kK covers the complete NLO terms. 
In Eq. (|D the <7 denote the leading order Wilson-coefficients, so that both the r t] and the 
kij are scheme-independent. PH is the tree-level phase-space factor given by 


(i-®i ^ds 

PH(ay, X2, X3) = 12 / —(s — x\ — rrolfl + x\ — s)w(s, xl, xl)w(s, x 2 , 1) (5) 

J s 

(X2+X3) 2 


with 

w(a , b, c ) = ( a 2 + b 2 + c 2 — 2 ab — 2 ac — 2 be) 1 / 2 . ( 6 ) 

The arguments of the phase-space factor are ratios of the quark masses, x c = m c /mb and 
x s = m s /mb • The weight functions f t] are tabulated in Table 1; they depend onQ 




/ = 


PH(x c 






w( 


)w(l,s,x 2 c ){s + x 2 -x 2 c )(l + s-x 2 c ). (7) 


( Xc+Xs ) 2 


/ describes the interference of operators with Dirac structure (V — A) 0 (V + A) with those 
of the form (V — A)®{V — A) and vanishes for zero final state quark masses. The coefficients 


r^ can be obtained from Ref. [|Tl|. In particular, we find 


2 i 


i= 1 j =1 


10863 - 1278n/ + 80n^ 

a s (mw) 

4//3 ° 15021 - 1530n/ + 80 n) 

a s {mw) 

162/12 

a s (/i) 

162/52 

Q^ s (/i) 


-8/0o 
. ( 8 ) 


3 Note that we have corrected a sign error in Ref. 
G. Buchalla for pointing out this mistake. 


|fOf in all the terms containing /. We are grateful to 
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where (3$ = 11 — 2n//3 is the lowest order coefficient of the QCD (3- function; in our case we 
have nf — 5 quark flavours, ku and k 22 were already given in Ref. |10|; k Vi can be obtained 
from the diagrams shown in Fig. 1 as 4 


k \2 = k 22 + - ( H e +B ) with H e PR(x c 
o 


7687T 5 

-^-^Im(VI+VIII+X+X t +XI+XI t ). (9) 
9s m b 


Here R is a scheme-dependent constant that removes the scheme-dependence of R e ; in naive 
dimensional regularization (see below) one finds B = 11 JO]]. Note that H e is independent 
of the definition of the quark mass, which only affects ku and k 22 through the self-energy 
diagrams. 

3. In the following, we present our results in the limit of vanishing strange quark mass 
(as will be discussed below, their dependence on this parameter is small); the full results 
are available from the authors as a Mathematica hie. We have checked that all formulae 
coincide with the corresponding ones in Ref. Q] when the appropriate limit is taken.[] 

Without going into too many details, we present first a short outline of the method of 
calculation of H e which is described at length in Ref. || . I 11 calculating the imaginary parts 
of the diagrams of Fig. 1, we use MS subtraction and control the ultraviolet divergences 
through dimensional regularization with an anticommuting 75 , often referred to as naive 
dimensional regularization (NDR). NDR is applicable if one uses Fierz-transformations to 
relate diagrams with closed fermion loops, which are ambiguous in NDR, to diagrams which 
are well-defined in NDR. As shown in Ref. j[T|| . Fierz-transformations are only valid diagram 
by diagram with a correct choice of the so-called evanescent operators. We have verified that 
in the limits m c , m s —> 0 our procedure yields the same results as obtained in other schemes 
0 . Technically, we calculate the imaginary parts of the forward-scattering amplitudes by 
applying Cutkosky rules. We regularize intermediate infra-red singularities by introducing 
small quark and gluon masses, denoted by p and A respectively, which allows phase-space 
integration to be done in four dimensions. 

For the sake of compactness in displaying the formulae, the square masses of the heavy 
quarks, c, b, are denoted by c, b. In the same spirit, we define A = (Vb — a/c) 2 ; w = w(b, c, t ) 
(or w = w(plpll 2 ) in Eqs. O. (B below), where w was defined in Eq. (^), and similarly 
v = w(c,c,t)/t. Finally, we omit the arguments of the functions. Our results are written 
in terms of the functions A, B , C, R, K, which were defined in Ref. ||, Eqs. (A.1)-(A.4), 
with M 2 = p\ and p 2 = p\. Hence, throughout this letter we consider the external square 
momenta p 2 , p\ and (pi +P 2) 2 to be the natural arguments of those functions. The integrals 
can be computed following standard techniques and be expressed in terms of logarithmic 
and dilogarithmic functions. The final analytic expressions for A, R, C, R, K are rather 
involved; we will give them elsewhere along with details of the calculation. The functions 
/C,-, j — 0 , 1 ,..., 7, denote certain phase-space integrals defined as 


(/Co,...,/C 7 ) = I LIPS(pi,p 2 , k) ^1, 2 p 1 k + \* ,Plk ' 


p 2 k 


p 2 k 


2 Pl k + X 2 ' A 2 -2 kV A 2 -2 kV 


4 We use the same notations as in Ref. [ 3|. 


5 Note a misprint in Eq. (C.9) in Ref. 


: the factor (2m 2 + s) should read (?n 2 + 2s). 
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(A 2 - 2kl)(2p 2 k + A 2 ) ’ (2pik + A 2 ) (2 p 2 k + A 2 ) 


, ( 10 ) 


where l = p\ + p 2 + k. /Co,- • ■, /C 5 , can be obtained from the results given in Ref. [14 


whereas the calculation of /Cg and /C 7 requires some effort. The prime symbol ('), as in C 
or K -2 below, will always denote the replacement p\ p 2 . After a tedious calculation, one 
finds 

r 2 


7T- 


(K«,W) = iI j(-C,[C + C']), 


( 11 ) 


where C is given by 

o, m iK+ + m 2 , K + l 2 

C = -2 In-In ——-— 

mi w(K + — 1 ) 


2 Ln 


' K + - 1 \ 

K + + m 2 /m 1 ) 


miK _ + m 2 . K + - AC 5 

— 2 In-in-— L 2 

mi 1 - AC 2 


2 In —- In K + — 2 L? 

mi 


l 2 — (mi + m 2 ) 2 

I 2 


-rri\K + \ 


+ 2A; 


AC-/4 \ 

AC + m 2 /mi J 

■mi 


-2A, 


AC - 1 


m 2 


m 2 


l 2 

+ In —- In 


A'- + m 2 /mi) 

l 2 


m 2 l 2 — (mi + m 2 y 


+ In 


AC+m 2 /mi A 


In ■ 


K + + m 2 /m 1 y/p 


— —Lr> 


' k + -k_ \ 

K+ + m 2 /mi J 


( 12 ) 


In the above formula, we have introduced the obvious notation y'p 2 = mj. The functions 
K± are given by 


AC = 


i 2 — pI~ pI±w 


(13) 


2mim 2 

Now that we have calculated the /Q phase-space integrals, we may give the imaginary parts 
of the relevant diagrams, where we denote the sum of all j-particle cuts for a given diagram 
by the superscript ^. To start with, we find for diagram VI 


Irn VI (j) = 


1 


87 lb J 4 c 

where the spectral densities p[ J> and p^P are given by 


d t (b - t) 2 \bp[ J - 2p^ 1} 


(14) 


M 2 ) =Re 4 

24? tH 

pP — Re j—j 
H2 48tt 4 


I (4c — t ) [t(A + B ) + 2(t + c)B] — 2(t + 2c) ^'C + — ^ + ( t 2 — 4c 2 )A , 
2 (t - c)[t(A + B) + 2C + 1 - (t - 2c)K] + (At 2 - 15 ct + 2c 2 )a1 , 


(3) 

Pi 


(3) 

P\ 


9 


67r 6 t 2 

J2 


9 


127r 6 t 


| t(t 2 - 4c 2 )/C 7 - 2 t(t + c)/C 1 - 2t/C 0 + 8 c/C 3 + 8/C 2 |, 

| - 2 t(t — c)(t — 2c)/C 7 + 2t(2t - c)Ki + tJC 0 - 4c/C 3 - 4/C 2 J . 


(15) 
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The external momenta in A, B , C, B, K satisfy p\ = p\ = c, (pi +P 2) 2 = whereas in p®, 
t = (jpi + P 2 + k) 2 . For the remaining diagrams, we obtain: 


Im VIII® = Re — z— / — (c - t) 2 (t - c - 6)tu 
327r 5 0 Jc t 


x <{ 2tA + 2TB - 2(c - t)B + 8 ( C + — ) + (c - t)K \, 


(16) 


Im VIII® = 


9 


2 ,A 


J dt (t — c — 5)u> |/C 0 — (t — c)/Ci — f/C', + (t — c) 2 /Cy j , (17) 


87 r 7 b 

g 2 /-A dt 


Im[x+xt](3) = 3 hbl 7 {t ~ c)2w r +c ~ 


2t(A + B) + 8C + 


+(b + c - t)k 


2 (t 2 — 2 tc — 2 tb + c 2 + k)B |, 


(18) 


Im 


[X + X*]« 4 > = jf y (t - c? {(* - C - bfKs + K 0 - (i - 6)K, - (t - c)X 4 } , 


(19) 


Im [XI + XI 1 ]® = —a— / — (6 - t) 2 n {(t + 2c)6 [t(A + IB) + {t + 2b) B 
1927 r b b J4c t 1 L 


- 2C - 1 - (6 - t)/C 

- 2t(t - c) [(t + 6)(A + B) - {b - t)(2B - K) + 2(7 + 1 
( 4 ) _ _&!_ f b 


}. 


( 20 ) 


lla[X1 + XlT = JL-j^ v { (t + 2c) 


tJCo + b(t — 6)/Ci + 2 /C 2 + b 2 } C 4 — 26/C 5 


- 6(4 - !)) 2 K 6 ] + 24(4 - c) [(f - 6)K, + t> C 4 - (t - 6) 2 K 6 ]} . 


( 21 ) 


In Eq. (U), the arguments of A, B , C, B, K satisfy p\ = c, pi = p 2 —> 0. As mentioned 
above, p (and also A) regularizes the infra-red singularities arising in intermediate steps of 
the calculation. They cancel upon addition of the 3- and 4-particle cut contributions to 
each diagram. In Eq. (|T8|), we set p\ — c, p\ — b. Finally, in Eq. (I2UJ), we have p\ — p 2 —> 0, 
p 2 = b, and in all three equations (pi + P 2) 2 = t. In Eq. (|T7| ), the arguments of the functions 
ICj are p\ = p 2 —> 0, p 2 = c and = t. In Eq. (|T9|), p\ = c, p 2 = t and l 2 = b. Finally, in 
Eq. (pl|), p\ = p 2 —> 0, p 2 = t and l 2 = b. 

The numerical results of our calculation are presented in Table 2, namely ku as a function 
of the charm quark mass for zero strange quark mass and m s = 0.2 GeV, respectively. For 
comparison and completeness we likewise give the coefficients k\\ and /C22 referring to the 
on-shell definition of quark masses. The table shows also the leading order correction k and 
the ratio Ynlo/Ylo — K. In the latter quantity, we have estimated the unknown NLO 
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penguin contributions very conservatively by assuming 0 < dij < 2 , which corresponds to 
\kij | < 15 for a s = 0.2. The other input parameters are given in the table caption. The effect 
of a finite value of the strange quark mass is tiny for the rate (though appreciable for the 
NLO corrections after dividing out the phase-space factor), and less than 5% for ay < 0.4, 
which is less than the estimated uncertainty from the unknown NLO penguin contributions. 
Using pole masses and a renormalization scale // = mj, we thus observe that T (6 —> ccs ) 
increases by (32 ± 7)% through NLO corrections for a reasonable choice of quark masses 
x c = 0.3. If we allow the renormalization scale to vary in the rangc(] m&/2 < /x < 2 mb and 
take the uncertainty in ay to be ±0.05, we obtain T nlo/ r lo = 1-32 ± 0.15. 

So far we have used the on-shell definition of the quark mass. However, far from being 
compulsory, this definition most likely introduces artificially large higher order perturbative 
corrections (cf. j4j). It is therefore most instructive to eliminate the pole mass in favour of 
an off-shell renormalized mass, such as, e.g. the MS mass. As discussed in Ref. [15 , this 
amounts to the replacement 


m b b PH (ay, x c , 0) —> m\ PH (ay, ay, 0) <! 1 ± 


a. 


7T 


20 


=,2 


— 5 In 


m 


A 


2 - 2 ay In x c 


dlnPH(x c , x c , 0) 
dx r 


( 22 ) 

in the decay rate, where x denotes a running quantity evaluated at the scale /i. We then 
obtainf] T nlo/ r lo = 1-2 ± 0.4, which indicates that the uncertainty due to unknown higher 
order corrections is appreciable. We shall come back to this point in the next section. 


4 . With the results for T(b —> ccs ) in hand, we are ready to give an updated value for 
the semileptonic branching ratio B$l of B mesons defined as 


Bsl = 


T(B —> Xev) 


E <= e )ftT r(B —> Xlvt) + T(B -»• X c ) ± T(B —> X c5 ) ± T(rare decays)' 


(23) 


Performing an expansion in the inverse b quark mass, it is possible to show |T| that the 
inclusive decay rate of a B meson into a final state X coincides with that of the underlying 
b quark decay up to corrections of order l/m% (n > 2 ): 


r (B -> X) = T(b -> x) (l ± 0(1/ml) ) . (24) 

The power-suppressed correction terms to the total inclusive widths of both semi- and 
nonleptonic decays were calculated in Refs. They depend on two hadronic matrix 

elements, Ai and A 2 . While the latter is related to the squared mass difference of the B and 
the B* meson, 

A 2 ~ ^ (m 2 B * — m 2 B ) = 0.12 GeV 2 , (25) 

6 We conform here to the conservative choice of “characteristic scales” that is preferred if one follows 
standard renormalization group improvement arguments, where large logarithms of type In m^//r 2 are to 
be avoided. The results obtained from summing the asymptotically leading part of the perturbative series 
seem, however, to indicate a lower scale, at least in semileptonic decays, cf. [|j, |j, [|]. It remains to be seen 
if those scales necessarily have to coincide or not. 

7 Note that x c = 0.30 ± 0.05 translates to x c (jj, = mb) = 0.28 ± 0.05. 
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the former, Ai, is difficult to measure; in this letter we use Ai = —(0.6 ± 0.1) GeV 2 as 
determined from QCD sum rules |[T7|. 

The one-loop corrections to the partial decay widths in Eq. (p3|) can be found tabulated 
in Ref. ]lOj except for T(fc —> ccs), which was a rough estimate. In Table |2] we give the 
corrections to this partial width. As in Ref. [II|, we neglect the rare decays in the present 
analysis because of their smallness. Using the same input parameters as in the last section, 
we find 


B sl = 12.0 ±0.9 


+0.9 

-1.35 


B sl = 11.2 ±1.0 


+1.0 

- 2 . 2 ' 


(26) 


Here the first error conies from the uncertainty in the input parameters x c , \i and T(6 —> ccs) 
(in which the effect of the penguin operators has been estimated), whereas the second one 
indicates the variation of the result with the renormalization scale (i. Both results are in 
agreement with the most recent experimental data and the particle data group world average 
B S l = (10.43 ±0.24)% [fT%| . Nevertheless, we observe a nonnegligible scheme-dependence of 
the two results. Although at the considered order in a s it is difficult to judge which scheme 
is “best”, we remark that at least for the semileptonic width one can sum up a certain 
class of terms, which are of order /3 qCk" + 1 . One observes that both, the explicit coefficients 
multiplying /3q ct” +1 with n not too big (say n < 5), and the resummed all-order expression 
are smaller in the MS than in the on-shell scheme (||], see also |^]). Interpreting this result 
with due caution, since the evidence that these terms are dominant already in low orders 
comes from empirical observation (of quantities with known complete ck 2 corrections) rather 
than from a theoretical principle, we still feel that it favours the MS scheme. Any further 
discussion would require the knowledge of complete ck 2 or even higher order terms, whose 
calculation is a formidable task. 

Finally, we would like to discuss shortly the average charm quark content of B decays, 
which is defined by 


(w) 


1 ± 


r(B 


Xcc) 


tot 


(27) 


We obtain (using again the same input parameters as in the last section): 


(n c ) — 1.27 ± 0.07, (fi c ) = 1.35 ±0.19, 


(28) 


which has to be compared with the experimental result ( n c ) exp = 1.04 ± 0.07 |I8|. The 
experimental and the theoretical numbers differ by 3 standard deviations. Unfortunately, we 
do not see any natural theoretical explanation for that fact, unless x c were much smaller than 
we assumed, which is however in conflict with the results obtained from the phenomenology 
of charmed particles. 
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Figure 



Figure 1: The diagrams contributing to the forward-scattering amplitude Eq. (Jl|) up to 
order a s without penguins. The crossed circles denote insertions of any of the operators Oi. 
Of the three internal quark lines, the upper one denotes the c quark, the lower one the s 
quark, and the middle one the c antiquark. 
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Tables 


fij 

1 2 3 4 5 6 

1 

1 

2 

2 1 

3 

2-1 

Z 3 1 

4 

2 2 2 i 

3 Z 3 1 

5 

2 / If 2/ §/ 1 

6 

If 2/ §/ 2/ | 1 


Table 1: Coefficients fij defined in Eq. (|3|). 


X c 

k(x c , 0 ,m b ) 


k 12 (n = m b ) 

^22 

K(x c , 0 , m b ) 

k K PH 

0 

1.054 

-1.34 

-7.75 

-1.41 

1.01 ±0.05 

1.065 ±0.059 

0.1 

1.052 

-0.14 

-6.31 

-0.53 

1.07 ±0.06 

0.959 ±0.052 

0.2 

1.047 

2.40 

-3.50 

0.99 

1.17 ± 0.06 

0.634 ±0.035 

0.3 

1.040 

6.44 

0.82 

2.99 

1.29 ±0.07 

0.263 ±0.015 

0.4 

1.032 

14.76 

9.50 

5.83 

1.45 ±0.08 

0.039 ±0.002 

x c 

hi(x c ,x s ,m b ) 

ku 

ki2 (h = m b ) 

^22 

K(x c ,x s ,m b ) 

k K PH 

0 

1.054 

-1.33 

-7.63 

-1.26 

1.02 ±0.05 

1.062 ±0.058 

0.1 

1.052 

-0.05 

-6.20 

-0.35 

1.08 ±0.06 

0.956 ±0.052 

0.2 

1.047 

2.53 

-3.36 

1.23 

1.18 ±0.06 

0.631 ±0.034 

0.3 

1.040 

6.69 

1.08 

3.41 

1.32 ±0.07 

0.259 ±0.014 

0.4 

1.032 

15.68 

10.43 

7.09 

1.54 ±0.08 

0.037 ±0.002 


Table 2: The LO and NLO corrections to the nonleptonic decay b —> ccs as a function of 
x c = m c /m b . The penultimate column gives the increase of the decay rate T(6 —*■ ccs ) in 
NLO if we include finite c and s quark effects in the radiative corrections. In the upper table 
we have put rn s = 0. The errors in K represent a conservative estimate of the unknown 
parts of the NLO penguin contributions. The input parameters are /i = mb = 4.8 GeV, 
x s = 0.04 and = 312 MeV, which corresponds to a s {mz) = 0.117. A comparison of the 
last column in both tables shows that the effect of the strange quark mass is negligible and 
actually below 5% for all values of x c . 
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